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Abstract 



We investigate the possibility that dark matter is made of heavy Dirac neutri- 
nos with mass m G [0(1) GeV- a few TeV] and with suppressed but non-zero cou- 
phng to the Standard Model Z as well as a coupling to an additional Z' gauge bo- 
son. The first part of this paper provides a model-independent analysis for the relic 
density and direct detection in terms of four main parameters: the mass, the cou- 
plings to the Z, to the Z' and to the Higgs. These WIMP candidates arise naturally 
as Kaluza-Klein states in extra-dimensional models with extended electroweak gauge 
group SU{2)l X SU{2)ji x U{1). They can be stable because of Kaluza-Klein parity or 
of other discrete symmetries related to baryon number for instance, or even, in the low 
mass and low coupling limits, just because of a phase-space-suppressed decay width. 
An interesting aspect of warped models is that the extra Z' typically couples only to 
the third generation, thus avoiding the usual experimental constraints. In the second 
part of the paper, we illustrate the situation in details in a warped GUT model. 
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Introduction 



It is well-known that a heavy (m > 1 GeV) Dirac Neutrino with Standard Model interactions 
is ruled out as dark matter because of its large coupling to the Z. On the one hand, it 
annihilates too strongly into Z to have the right thermal abundance. On the other hand, 
even if it had the right relic density from a non-standard production mechanism, it would 
scatter elastically off nuclei with a large cross section induced by the Z exchange and should 
have been seen in direct detection experiments, unless its mass is larger than several tens of 
TeV [H [2] . Moreover, this type of neutrino is excluded by electroweak precision tests [3] . In 
contrast, a sterile Majorana neutrino as dark matter is a possibility that has raised interest 
lately [HE]. In this case, the neutrino mass that has been considered is rather in the keV- 
MeV range and behaves as warm dark matter if < 20 keV. In addition, the Majorana 
mass scale is determined by the see-saw formula for neutrino masses. 

In the present work, we are considering a different type of Dirac neutrino, denoted v\ 
corresponding to a typical cold dark matter WIMP candidate with a mass at the electroweak 
scale, suggesting that dark matter and the electroweak scale are somehow related (even if 
the mass of the dark matter particle does not come from electroweak symmetry breaking). 
We do not assume any particular relation between the mass scale of v' and that of the light 
standard model neutrinos. The relic density of v' is entirely determined by the standard 
thermal mechanism and therefore by its annihilation cross section. We are assuming that 
the coupling to the Standard Model Z is suppressed. There are various reasons why this 
can happen. A typical framework is to start with an SU{2)l singlet neutrino but charged 
under SU{2)r. Because the gauge bosons of SU{2)r are heavy, their interactions with are 
quite feeble and this makes vr behave as a WIMP. In addition, electroweak (EW) symmetry 
breaking typically induces a mixing between Z and Z' , leading to an effective small coupling 
of vr to the Z. Examples of this type were studied in warped extra dimensions [SI El E], 
and in universal extra dimensions [S]. Note that the most promising realistic warped extra- 
dimensional scenarios need the EW gauge group to be extended to SU{2)l x SU{2)r x U{1). 
In this context, Kaluza-Klein Dirac neutrinos charged under SU{2)r are necessarily part of 
these constructions, even though their stability typically requires an additional ingredient. 
For instance, it was shown in [HI [7| that implementing baryon number conservation in warped 
GUTs leads to the stability of a KK RH neutrino. 

Finally, even in the absence of an additional symmetry, z/' can be cosmologically stable if 
the couplings involved in its decay are very suppressed. This can happen even if the neutrino 
has a large annihilation cross section providing the correct relic density. 

Dirac neutrino dark matter was also studied in the 4D models of Ref. [TU] where the 
Z-coupling suppression has a different origin and results from mass mixing between gauge 
eigen states with opposite isospin. Ref. [TDJ did not consider the effect of a Z' that we study 
in details here. In the first part of this paper, we present a model-independent analysis for 
the viability of Dirac Neutrino dark matter in terms of three main parameters, the mass, 
the Z coupling and the Z' coupling of u'. We also discuss the effect of a coupling with the 
Higgs. The remaining part of the paper is a refined analysis of the dark matter candidate 
which arises in the warped GUT models of Ref. P, In both studies, the computation of 
the relic density is performed with micrOMEGAs2.0 [Tl] after implementing new model files 
into CalcHEP [12]. 
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Part I - Model-independent analysis 



In this part, we consider a generic extension of the Standard Model (SM) containing 
a stable heavy Dirac neutrino denoted u', with mass M^/. We assume a parity symmetry 
(except in Section \2.2\i that does not act on SM fields and u' is the lightest new particle 
charged under it. The model also contains an additional Z' gauge boson and potentially 
a charged gauge boson W, with masses Mz', > 500 GeV. To avoid most low energy 

constraints we will assume that these new gauge bosons couple only to the fermions of the 
third generation. 6zz' is the mixing which induces the Z'WW and Z'ZH couplings. We 
also introduce gn as the Higgs coupling to the z/' (for instance induced via mixing with a 
heavy as illustrated in Part II). We assume that only one chirality of u' couples to the 
gauge bosons (in our numerical examples, we chose the right-handed chirality). The effective 
couplings of u' to Z, Z' and H are denoted gz, gz' and gn respectively: 



We work at the level of a low energy effective theory. We assume that the remaining new 
physics which makes the model more complete does not interfere much with our dark matter 
analysis. 

1 Direct detection constraints 

Direct detection constraints are very simple for a Dirac neutrino (see for instance section 
3 of Ref. |1]). The cross-section ajsf for u' scattering on nucleons is governed by the t- 
channel exchange of the Z, the Z' exchange is comparatively negligible. In contrast with 
Majorana dark matter, the Z exchange contributes to the spin-independent scattering cross 
section. Therefore, strong constraints on gz, the u' coupling to the Z, are derived from 
direct detection experiments, in particular CDMS |13] and recently XENON which has now 
the most stringent limit [Hj. The dependence on M^/ that we observe for this constraint is 
related to the experimental sensitivity, which is optimal around 50-100 GeV. The theoretical 
prediction does not depend on M^/, as illustrated on Fig.dl In our plots, we use the parameter 
g/gz where g is the SM electroweak coupling g = e/{sm9w)- The direct detection limit has 
been rescaled to take into account the fact that, for a Dirac fermion, the interaction of the 
Z with protons is suppressed by a factor (1 — 4 sin^ 9wY (see Eq. |2]) so that the scattering 
with nucleons is completely dominated by neutrons. When CDMS and XENON quote their 
limit they rather assume that protons and neutrons contribute equally. This means that the 
CDMS and XENON exclusion curves go up by a factor AV((1 -Asin^ ew)Z - (A - Z)f 
i.e 3.7 and 3.4 respectively. We should keep in mind that the bound from direct detection 
experiments is subject to some astrophysical uncertainties such as the velocity distribution 
of the WIMP. Based on Ref. [15] we could estimate these uncertainties and allow for a factor 
~ 3 in the interpretation of the CDMS and XENON limits, over the full 10 GeV - 1 TeV 
WIMP mass range even though there are actually much larger uncertainties for masses below 
~ 40 GeV. For clarity, we have not displayed this uncertainty in our plots. 

It is clear from Fig. [T]that we have to impose gz ^ fi'/lOO for M^i < 400 GeV to satisfy 
the XENON constraint. If v' has a sizable coupling to the Higgs, the elastic scattering 
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Figure 1: Neutrino- neutron scattering cross section due to Z-exchange for q/qz = 
10,30,100,300 where g = e/(sin6'vK) is the SM couphng. The dotted hne shows the effect 
of adding the Higgs exchange for gH=0.25, rriH = 120 GeV, in the case where g/g^ = 300. 
Also represented is the CDMS hmit as well as the recent XENON limit [H]. 



via Higgs exchange is not always negligible compared to the Z-exchange especially when u' 
couples weakly to the Z. The spin-independent elastic scattering cross section on nucleons is 
the sum of two contributions (when averaging over z/' and z/', the negative interference term 
cancels) : 
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where /^'^ are taken from Ref. [16]. For example, when mu = 120 GeV, gH = 0.25 and 
g/ gz = 300, the Higgs and Z contributions become comparable for elastic scattering on neu- 
trons, as illustrated in Fig. [H while for elastic scattering on protons, the Higgs contribution 
dominates. 

In contrast with Majorana fermions like in the MSSM, the parameter that determines 
the elastic scattering of Dirac neutrino dark matter on nucleus is the same parameter that 
drives annihilation and determines the relic density. We now look at when Dirac neutrinos 
can inherit the correct thermal abundance. 
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2 Annihilation 
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Figure 2: v' annihilation diagrams into //, WW, Zh, hh and ZZ. 



2.1 Annihilation via s-channel Z exchange 

Annihilation channels of z/' are listed in Fig. O We first look at the effect of the coupling to 
the Z, gz, and set gz' = 0. For M^r < 100 GeV, annihilation into fermions dominate due 
to the Z resonance. For M^/ > 100 GeV, main annihilation channels are into WW and Zh 
(the relative contribution of Zh increases for larger M^i). The contribution from top pairs 
is small. Figure [3^. shows how VLyiK^ decreases as function of Myi for different gz couplings. 
In Fig. [3^, we show for comparison the prediction for a fourth generation Dirac neutrino 
with SM coupling to the Z and a Yukawa coupling to the Higgs as studied in [IT]. Since v' 
has non-standard couplings, the total annihilation cross section grows with MJ^ Unitarity 
breaks down for M^i at the multi TeV scale and we do not show any predictions beyond 
these values. Figure [3)d shows the effect of a coupling of v' with the Higgtl, in particular the 
resonance at M^i ~ Mh/2. The contributions of the different annihilation channels are also 
displayed in Fig. HI The ZZ channel is important if v' has a sizable coupling to the Higgs. 

Near the Z resonance, only a weak coupling is necessary to get VL^ih^ ^ 0.1. We combine 
the relic density constraints with the direct detection constraints in Figure E] where we show 
the WMAP [in] allowed region (0.097 < Vth^ < 0.113)1 in the g/gz - M^, plane, as well as 
the CDMS and Xenon limits. The region satisfying both constraints corresponds to either 
M^/ ^ 40-50 GeV or M^/ > 500 GeV if the Higgs coupling is negligible. For M^, > 500 
GeV, the relic density constraint on is more severe than the direct detection constraint. 
If the Higgs coupling to u' is sizable, a region around M^/ ~ Mh/2 opens up. 



^This is different from the behaviour of the vector- Hke Kaluza-Klein neutrino studied in [TH]. 

^We neglect the contribution to the mass of i^' due to the Yukawa couphng. 

•^A more conservative upper bound, 0.097 < flh^ <0.134, was recently advocated in [SUj . 
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Figure 3: n^'h^ versus u' for g/g, = 10,30, 100,300, 1000 where g = e/isinOw) is the SM 
couphng. In a), only the Z-exchange is included. For comparison, we also show in dotted 
line the relic density of a fourth generation Dirac neutrino with SM coupling to the Z and 
a Yukawa coupling to the Higgs. For this curve, we set rriH = 200 GeV. In b) the Higgs 
exchange is included as well with gH=0.25. 




Figure 4: Annihilation cross sections at freese-out. The effect of a Z' is omitted here and 
will be shown in Fig. 
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Figure 5: WMAP region, 0.097 < fi/i^ < 0.113, in the {Myi,g/gz) plane (where g = e/sinOw 
is the SM couphng). The blue band includes the Z-exchange only, the red band includes the 
Higgs exchange as well. The region excluded by CDMS and Xenon is also displayed. The 
dotted line takes into account the Higgs contribution. The region above the direct detection 
line and below the WMAP band is allowed. In this figure, the effect of the Z' is omitted. 



2.2 Light neutrino scenario 

For WIMP masses below 10 GeV, direct detection constraints do not apply. We concentrate 
on this particular case in this section. Light WIMPs also potentially offer the interesting 
prospect that they can be stable without the need to introduce an extra symmetry, just 
because there is no SM particle they can decay into, at least on a cosmological time scale. 
In the absence of any particular discrete symmetry, our exotic light neutrino can have a 
three-body decay, like the muon. For instance, if the underlying theory is 5-dimensional 
with SU{2)ji gauge symmetry, there are gauge couplings between the KK mode and the 
zero mode of ur: v^^^Z'vf, v^r^WrI^. They induce, via Z - Z' and Wr-W mixings, the 
effective couplings {v' then corresponds to v'^r)'- Zu'iy^ and Wv'l^. If these are sufficiently 
suppressed, with strength of order g^i, the lifetime of v' can indeed exceed the age of the 
universe, as illustrated in Fig. [6l using the width F ~ {g^, / gfG'j.Ml, / {192ti^). For a neutrino 
mass in the 1-10 GeV range, these coupling should be smaller than 10~^^ for the neutrino 
to be cosmologically stable. This is actually naturally realized in Randall- Sundrum models 
where the zero mode neutrino is localized near the Planck brane while the KK vr is peaked 
on the TeV brane. The overlap of their wave functions is therefore very suppressed, resulting 
in a tiny effective 4D coupling. However, in this case, we expect the KK vr to have a mass 
in the TeV range rather than 10 GeV. The only KK fermions that can be naturally light 
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are those belonging to the multiplet containing the top quark, as discussed in Part II. For 
instance, in warped GUT models, the KK belonging to the 10 of S'O(IO) which contains 
the top quark (that is the only zero-mode SM particle in the 10) does not couple directly 
to any light SM fermion. The only SM fermion it can directly couple to is the top, via 
an 5*0(10) gauge boson. Therefore, its decay has to go through a very large number of 
intermediate states and will be very suppressed. In other words, the boundary conditions on 
the different components of the multiplet containing are such that no three-body decay 
is allowed. 

Now the question is whether such a neutrino can naturally inherit the correct abundance. 
The interesting aspect here is that these models typically offer the possibility that the multi- 
body decay is suppressed while the self-annihilation can be large. This can be explained in 
terms of the different localizations for the wave functions of SM light fermions on one side 
and KK modes (of both fermions and gauge bosons) on the other side. 

M„ (GeV) 




Figure 6: g^i is the typical size for the couplings Zv'v and Wv'l^. Below the line, is 
cosmologically stable. 

It is clear from Fig. [3^ that for M^i ~ O(GeV), gz ~ 0{1) is needed to obtain the 
correct thermal relic density; this is in contradiction with the experimental constraint from 
the invisible decay width of the Z which requires gz ^ 0.035. One way to open the window 
could be to consider that v' couples to a light ((9(GeV)) singlet scalar field (p that decays 
into SM fermions via its mixing with the Higgs. Near the resonance, irtyi ~ the right 

amount of annihilation can be obtained. Another possibility is to assume that the reheat 
temperature is below the freese-out temperature Tp ^ M^'/25. For M^/ ~ 10 GeV, Tp ~ 0.4 
GeV, this is still well above the BBN temperature. In this case, the neutrino is produced 
through scattering in the plasma and the correct relic abundance may be reproduced with 
an appropriate choice of reheat temperature and gz coupling. Note that in addition to u' 
production via the gz, gz' and gn couplings, there is production of via the g^' couplings 
defined above, which is similar to production via neutrino oscillation [H |2I] but it is very 
suppressed given the g^' values that we consider for cosmological stability. 

Even if the lifetime of the neutrino exceeds the age of the universe, there are additional 
bounds on gi,/ from EGRET and COMPTEL measurements of the diffuse gamma ray spec- 
trum, that constrain the radiative decay with gamma emission u' — * wy [22]. In conclusion, 
this 0{GeV) Dirac neutrino scenario would deserve a detailed analysis. 
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2.3 Annihilation via t-channel heavy lepton exchange 



We now consider the effect of the t-channel exchange of a heavy charged lepton r' (potentially 
a Kaluza-Klein lepton, see also the toy model of Ref. [10]) leading to the annihilation into 
W. The coupling is defined as 



where gw is written here in terms of the W — W mixing resulting from EW symmetry 
breaking. According to Fig. [TJ quite a large gw coupling is needed to see the effect of these 
diagrams. For instance, if gw arises from W — W mixing and Mw' = 1 TeV, then we need 
Xw ^ 10. We can make the following observations: 

• There is a destructive interference between the s and t channels leading to an increase 
in Vlyih'^ compared to the case with Z exchange only. 

• When one W can be produced in association with a in the annihilation, there is a 
decrease in Q^rh'^ (see the kink around M,^/ ~500 GeV). 

• When t' is not too much heavier than z/', coannihilation effects can reduce the relic 
density. For illustration, we show the case M^//Mj,/ = 1.2 where coannihilations more 
than compensate for the increase of Vt^ih'^ due to the t-channel exchange. In addition, 
near the W resonance, there is a sharp drop in the relic density. 



Figure 7: Effect of the t-channel annihilation on Qi.rh'^ for g/gz = 100. The blue full curve 
includes the Z s-exchange only, the dashed (dot) curve includes the t-channel with Mw' = 
1 TeV, Xw = 10, leading to gw = 0.065 and Mr' /Myi = 1.2 without (with) coannihilations. 




(3) 




M„(GeV) 
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2.4 Constraints on the W 



If the W comes from a Left-Right model, there are strong hmits on the W mass and mixing 
assuming a manifest LR symmetry i.e. the same mixing matrix in L and R quark sectors. 
Typical limits on the mixing angle are approximately 6 < 10~^ (assuming = g), while 
limits on the mass are around 1 TeV. The best limit is from — Ks mixing, M\yi > 
1.6 TeV [23]. There are also Tevatron limits using the decay channel into electron and right- 
handed neutrino leading to Mw' ^ 750 GeV [21]. In our analysis, we consider that only 
the coupling of W to the third generation is non-suppressed, as well-motivated in Randall- 
Sundrum models. In this case, the Tevatron constraints are weakened if we restrict the 
decay to the r channel and indirect limits from K or decays do not apply. Nevertheless, 
if W couples to quarks of the third generation (as in Randall-Sundrum models) there is 
an important constraint from b ^ sj [251 [26] . The W leads to an enhancement of the 
b —>■ s'-f amplitude by a factor mt/rrib- Irrespective of the W mass, a limit on the mixing is 
—0.015 < 6 X Qw'/g < 0.003 [25l 126] . This assumes similar quark mixing matrices in the L 
and R sectors. If the W coupling to leptons is similar to that of quarks, the above constraint 
is incompatible with the values chosen in Fig. [TJ xw = 10, M^/ = 1 TeV, corresponding to 
gw' = 0.065. Relaxing the assumption that the mixing matrix is the same in the left and 
right sectors will not help sufficiently. In conclusion, the W effects can be ignored. 



2.5 Annihilation via s-channel Z and Z' exchange 

We now consider the combined effect of the two annihilation channels through Z and Z'. 
We first look at the Z' coupling to W pairs, arising from Z-Z' mixing: 

gz'ww = 9zww ^zz' (4) 

where gzww is the Standard Model coupling. Next we add the interaction of the Z' with top 
quarks. There can be a noticeable reduction of the relic density due to this extra channel. 
Figure [HI shows the effect of the Z' exchange on ^yiln^ for different values of the parameters. 
It opens a region at the Z' resonance, for M^' ~ Mz' /2, and also reduces the upper bound 
on Myi from WMAP compared to the case with Z exchange only. We show in Fig. [9] the 
contributions of the different channels to the total annihilation cross section. There is a 
destructive interference between the Z and Z' contributions in the WW channel beyond the 
Z' resonance. However, this does not produce a significant increase in the relic density, since 
the Z' annihilation into Zh is still important. In Fig. [9] we chose gzz'H = 0.5gz'V. 



2.6 Constraints on the Z' 

The coupling of u' to the Z can be induced, for example, via Z-Z' mixing or via mixing with 
another heavy neutrino which has a large coupling to the Z. If the Z-Z' mixing comes from 
the Higgs vev only, then 

ozz' - - — ^ 9z' jjr 

2 cos Uw M^, 

In the absence of a custodial symmetry protecting the p (or T) parameter, we have to impose 
the constraint aT ^ 6'^^,M'^,/M'^ < 10~^. When the gz coupling is induced only via Z — Z' 
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Figure 8: Q^'h"^ for g/gz = 100, gu = ^ and gz' = (black), gz' = 0.3 (red). Also shown 
is the case where Z' couples in addition to all third generation fermions (green dash) with 
gz'ff = gz' = 0.3 in a) and gz'ff = 1 in b). We also show in b) the case with gz' = 1 (blue) 
and the case where Z' couples to top quarks with gz'tt = 1 (pink/dash-dot). Note that we 
fixed 9zz' = 10~^ in both a) and b) even though it is expected to scale as M|/M|,. 




Figure 9: Annihilation cross-sections at freeze-out into WW (solid), Zh (dash), ZZ (blue 
dash-dot), hh (red dash-dot) and ti (dot) for gz = g/lOO, gn = 0.25, niH = 200 GeV. In 
black, only the Z exchange is included {gz' = 0), in green, the Z' exchange is included with 
fl'z' = 1, gz'ww = ^^~^gzww and gzz'H = 0-5gz'V. 
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mixing, then Qzz' = Qz/qz' and it is difficult to satisfy the constraint on the T parameter 
if g^i ~ gz'- A similar LEP constraint comes from the shift in the vectorial coupling of the 
Z to the r and the 6, leading to 6zz' Qz'ff < 4 10~^. The T parameter constraint can be 
relaxed in models with SU{2)b. gauge symmetry which are in any case a strong motivation 
for Dirac neutrino dark matter. The constraint on the Zhh coupling can also be evaded (see 

e.g. m)- 

Like the limits on W , most direct collider searches involve fermions of the first and 
second generations and if we assume that Z' couples to the third generation only, we can 
tolerate Z' as light as ~ 500 GeV A which has generation-dependent couplings (like in 
Randall-Sundrum models) will induce tree-level FCNC. If it couples to the third generation 
only: 

>C = {gz'bLhf^Pib + gz'bRht^PRb)Z' (6) 

flavour non-diagonal couplings to the down-type quarks are induced 

^ = {gz'hJilp^PLUl^iUl^jdj + gz'bndi'lf^PRU*R3iUR3jdj)Z' (7) 

with Ul,r the mixing matrix for the left-(right)-handed down-type quarks. This is turn will 
induce a flavour non-diagonal coupling to the Z due to Z — Z' mixing. FCNC effects due 
to a Z' with non- universal couplings are analysed in Ref . [2HI [22] • Constraints are model- 
dependent and can be avoided. They are typically weaker than the T-parameter constraint. 




Figure 10: WMAP region in the [Myi , g / gz) plane with Z exchange only (blue), with Z' 
echange for gz' = 0.3 (green) and gz' = 1-2 (red). The effect of adding a Z'it coupling is 
also shown (pink) assuming g§,^^ = gz' = 1.2. In these plots, 9zz' = 10~^ and gn = 0. The 
allowed region is above the XENON dashed line and below the Qh"^ bands. 
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2.7 Contour plots 



We summarize our results and show the WMAP region in three different planes: {Myi , g / gz) 
in Fig. m and [IH {M^'.gz') in Fig. [HI and {M^,,Mz>) in Fig. [H In Fig. [12] and [131 
we neglect for simplicity the coupling of Z' to the fermions. In Fig. [TU] and [TT], we fixed 
Gzz' = 10~^ independently of the relation ^ while the two plots of Fig. [T2] satisfy Eq. 
In addition, the two right-handed plots of Fig. [12] assume that the only source of the gz 
coupling is the Ozz' mixing. As a result, the v' — Z coupling is suppressed and falls 
within the WMAP range only for a large coupling to Z\ gz' ~ 1, or for M^i M^//2. 
The situation is best summarized in Figure [11] which captures what is the allowed region of 
parameter space after imposing the WMAP bound and XENON constraints: There is a small 
mass window allowed for M^i ^ Mz/2. There are other wider mass windows near Mh/2 
and Mz'/2. Away from these resonance effects, a large region opens up for M^/ > 700 GeV. 




Figure 11: Summary plot: WMAP region in the {M^',g/gz) plane with Z exchange only 
(blue), with Z' echange for gz' = 0.3 and Mz' = 1 TeV (green) and Higgs exchange with 
gn = 0.5 (red). In this plot, Ozz' = 10"^ The CDMS and XENON direct detection 
constraints are shown with (dashed) and without (solid) a z/'-Higgs coupling. The allowed 
region of parameter space is above the XENON line and below the Qh^ bands. 
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Figure 12: WMAP region in the {Mt,',gz') plane. In all cases, 6zz' = 0.365'f,('y/M^/)^ and 
9z' — 9z'- In the two plots on the left, gz is fixed to gz = g/^OO (red) and gz = (7/300 (blue) 
independently of the gz' value. In the two plots on the right, gz is fully due to Z-Z' mixing. 

3 Collider signatures 

Like in other WIMP models, the standard searches rely on pair production of the heavier 
exotic particles which ultimately decay into the WIMP, leading to signals with energetic 
leptons and/or jets and missing Et- We list below some signatures which are more specific 
to the neutrino WIMP model. 
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Figure 13: Relic density in the {Mi,',Mz') plane for 9zz' = 10"'^, gn = dz' = 0.7, mn = 
115 GeV and we neglect the coupling of the Z' to fermions. The blue (dark grey band) 
corresponds to the WMAP range while in the green (light grey band) Q^ih'^ > 0.113 and in 
the white region Qiy'h"^ < 0.097. On the first plot, the contours of ^l^'h"^ = 0.097,0.113 for 
9z = fi'/lOO and 9zz' = 9z/ 9z' (dash) are also shown. 

3.1 Invisible Higgs decay into u' 

As seen previously, if v' has a significant coupling qh to the Higgs, it can account for the dark 
matter of the universe for m,^/ ~ mnj'l corresponding to the Higgs resonance. A significant 
gu coupling can arise for instance in the model of [HI Ej (see section 9.3). 

As a result, the Higgs can decay invisibly into v' with a significant branching fraction. 
The gu coupling could be probed at the LHC in Higgs production associated with gauge 
bosons [30], or with top quarks or in the weak boson fusion process [31|. The weak boson 
fusion process seems to be the most promising. In the following we use the results of Ref. [32j 
where an analysis including a detector simulation was performed and a limit on the invisible 
width was obtained for various Higgs masses. In general, this limit takes into account the 
fact that the production cross-section could be modified relative to the SM one, here we 
assume that the couplings of the Higgs to quarks are the standard ones. A more recent 
study [33] performed a refined analysis of the ZH channel, combined with the boson fusion 
channel. For light Higgses {vnn < 160GeV) the results are similar to the ones of Ref. 
The partial width of the Higgs into neutrinos is 

r(/.^.V) = f^M^(i-4^ 
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Fig. [m shows the resulting hmit that can be obtained on the qh couphng at LHC 




Figure 14: Above the hues, an invisible Higgs would be probed at LHC with C = 100 fb 



mtau (GeV ) 




Figure 15: Lines delimiting the region where r' is too long-lived to decay inside the detector 
for three values of A, the relative mass splitting between u' and r'. 



3.2 r' production 

Let us consider models where u' belongs to a gauge multiplet, for instance an SU(2)ji mul- 
tiplet. If u' and its r' partner are nearly degenerate in mass, r' is long-lived and the pair 
production of r' could lead to interesting stable CHAMF0-like signatures. This is to be 
contrasted with the standard scenario in which the dark matter is dominantly produced 
through decays of colored particle and therefore accompanied by energetic jets rather than 

''Charged Massive Particles 
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Figure 16: pp t'^t' cross section at the Tevatron and the LHC. M^/ = 3 TeV (sohd 
hne) and Mz' = 800 GeV (dashed hne). The couphng of Z' to b is gz'b = 1- 

by charged tracks. This situation was also discussed in the dark matter model of |34j . Fig. [TSl 
shows the region in the (fi'vyrV; "^r') plane where r' decays outside the detector. Limits from 
LEP [33] (m,-/ > 100 GeV) and DO [SB] are reproduced on our plot (Fig. [12]) of the r' pair 
production cross section, which is dominated by the Z exchange since Z' only couples to 
the third generation. A more likely possibility is that r' will decay inside the detector. In 
this case, the search is similar to that of sleptons with signature two leptons +E!^^^^ (see 
Ref. [371 EH] for LEP constraints). 

4 Conclusion of Part I 

In summary, Dirac neutrinos are viable dark matter candidates. The situation is best sum- 
marized in Figure [TH For a mass between 10 GeV and 500 GeV, the main requirement is 
that the coupling to the Z should be at least 100 times smaller than the SM neutrino-Z 
coupling in order to satisfy the direct detection constraint. Once this is satisfied, there is 
a large range of neutrino and Z' masses as well as Z' couplings that lead to the correct 
thermal abundance. The annihilation via Z is the dominant mechanism for z/' masses below 
100 GeV. Near Mz/2, the annihilation mechanism is even too efficient. If u' has a large 
coupling to the Higgs, M^i ~ M^/^ can lead to the correct relic density. Finally, u' — Z' 
couplings open a large spectrum of M^r possibilities in the multi hundred GeV range up to 
the Z' mass. 

In this work, we have assumed that there is no primordial leptonic asymmetry in the 
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dark matter sector {n^' =n^'). Obviously, our predictions for the relic density could change 
significantly if there was such an asymmetry, like there is in the visible matter sector. In 
constrast, this issue does not arise with neutralino dark matter or heavy Kaluza-Klein gauge 
boson dark matter. We have also restricted our analysis to the case where Dirac neutrinos 
would constitute all the dark matter. Constraints would be relaxed if we assumed instead 
that they constitute only a subdominant piece of dark matteiE]. Moreover, we only studied 
the case where Dirac neutrinos are thermal relics. Very different conclusions can be drawn 
if instead, the production mechanism is non thermal, and this is left for a future project. 

Except in Section 12.31 we have not considered coannihilation effects, as this is a more 
model-dependent issue. They will be studied in the explicit example of Part II. Finally, we 
have not discussed indirect searches in this work, this was done in [S] for the case of the LZP 
model that we now present in details. 



Part II - An explicit example: The 
LZP in warped GUTs 




In Ref. [HI [7] it was shown that in models of warped extra-dimensions embedded in a GUT, 
the symmetry introduced to prevent rapid proton decay, a Z3 symmetry, also guarantees the 
stability of a light KK fermion, a KK right-handed neutrino. This particle is called the LZP 
and its properties have been studied in [61 [7]. A detailed analysis of the indirect detection 
prospects in neutrino telescopes, cosmic positron experiments and gamma ray telescopes was 
also presented in and in [lO] for antiproton experiments. Some collider signatures were 
discussed in [7] and more recently in [H]. In this paper, we revisit the properties of the LZP 
and perform a complete calculation of its relic density. 

The underlying model is based on the Randall-Sundrum setup [42j, where the hierarchy 
between the electroweak (EW) and the Planck scales arises from a warped higher dimen- 
sional spacetime. All Standard Model (SM) fields except the Higgs (to solve the hierarchy 
problem, it is sufficient that just the Higgs -or alternative dynamics responsible for elec- 
troweak symmetry breaking- be localized at the TeV brane) have been promoted to bulk 
fields rather than brane fields. EW precision constraints require the EW gauge symmetry 
in the 5-dimensional bulk to be enlarged to SU{2)l x SU{2)r x U{1)x [IS]. The AdS/CFT 
correspondence suggests that this model is dual to a strongly coupled CFT Higgs sector [44j . 
Also, the SU{2)l x SU{2)r gauge symmetry in the RS bulk implies the presence of a global 
custodial isospin symmetry of the CFT Higgs sector, thus protecting EW observables from 
excessive new contributions [13] . 

In this framework, Kaluza-Klein (KK) excitations of gauge bosons of mass M^k ~ 
3 TeV are allowed and interestingly, light KK fermions are expected in the spectrum as a 
consequence of the heaviness of the top quark. The heaviness of the top quark is explained 
by the localization of the wave function of the top quark zero mode near the TeV brane. 
This is done by choosing a small 5D bulk mass (the so-called "c" parameter), guaranteeing 
a large Yukawa coupling with the Higgs. It is clear that if we want to embed these models 

^The scenario where a fourth generation neutrino with mass near the Mz/'2. window is a subdominant 
component of dark matter was studied in j39j . 
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into a GUT, we cannot consider SU{5) but rather Pati-Salam or 5*0(10) schemes. When 
the Right-Handed (RH) top quark is included in a GUT multiplet, its KK partners do not 
have a zero mode but their first KK excitation turn out to be hght. The masses and some 
of the couplings of these KK fermions are determined by the c-parameters (see Eq. [22]in the 
appendix). To have 0{1) top Yukawa, the right-handed top must have < 0. We fix the 
parameters associated with the top quark to be = cl = 0.4, = c/j = —0.5. Among all 
SM particles, the c of the RH top quark is the smallest. As a consequence, the KK modes 
inside its multiplet are predicted to be light (see Eq. [22]). They are likely to be the lightest 
KK states in these models and could be produced at the LHC [3 HI] . 

The main feature of unification in extra dimensions is that the SM fermions have to 
be split into different GUT multiplets (unless the SM fermions are localized on the Planck 
brane). However, this is still not enough to prevent proton decay from dangerous higher 
dimensional operators localized on the TeV brane. In P, [7] , the problem of baryon number 
violation was solved by imposing a symmetry. The consequence of the symmetry is 
also to provide a stable particle. In this model, for each fermion generation, there are at 
least three 16 multiplets. Each of them is assigned a baryon number. For instance, for the 
third generation: 

B= 1/3 ti^,hi^,t%,t%,UL,rL,T'^,UR 
B= -1/3 ti, 6l, t^, b^, i/L, tl, 4, 4 

B= tL,6L,t^,6^,z/L,rL,r^,z/R (9) 

where the particles in bold have zero modes and correspond to the SM fermions. The Z^ 
symmetry is 

$ ^ e^<^-^)^ (10) 

where ric (ng) is the number of color indices of the (anti)particle $. SM particles do not 
carry any Z^ charge. The Lightest Z3 charged Particle (LZP) is therefore stable. 

Although there are many new fermions in this model, we focus on the lightest ones, i.e. 
the level one KK fermions that belong to the multiplet containing the SM tji (the B = —1/3 
multiplet of the third generation in Eq. The LZP belongs to this multiplet. The only 
phenomenologically acceptable choice is that is the LZP (i/£, couples too strongly to the 
Z). All fermions inside a given multiplet should have the same c-parameter. However, 
because of bulk GUT breaking effects, there can be effectively large splittings so that the 
c-parameters of each component in the multiplet can be treated as free parameters. The 
c-parameters of members of SU{2) doublets are identical, for example Ct^ for both t^l\b^l\ 
The model therefore includes six free parameters for the third generation KK fermions. 

We ignore the KK partner of the right-handed top. has different boundary conditions to 
provide the zero mode of the RH top quark and its first KK mode is in the multi-TeV range. 

We also consider the influence of the of the i? = 1/3 multiplet which plays a role 
because of the mixing induced to the LZP. We will not include this particle in the model or 
compute explicitly its contribution to scattering processes but we will take into account its 
mixing with the LZP which will influence the Z-LZP coupling. 
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Among the 45 gauge bosons of 5*0 (10), we only consider the gauge bosons of SU{2)i x 
SU{2)r X U{1), W^, Wl, W^, W^, B' which neutral components recombine to give Z, Z' , 7, 
and a leptoquark gauge boson of electric charge Q = 2/3 which belongs to a color triplet + 
anti-triplet invariant under S't/(2)'s, XsiXs. This state has Z3 number ±1/3. Other gauge 
bosons will not directly enter the 2 — 2 annihilation cross-sections that are relevant for us. 
The masses of these KK gauge bosons are taken to be all equal, = = = Mkk- 

The parameters that we expect to have an effect on the annihilation rates are: ^'lo , the 
4D S'O(IO) coupling which determines the strength of the Z' couplings (thus the coupling 
of the Z to the LZP) and can be considered as a free parameter; Ci^^ which determines the 
LZP mass; Cp^ which enters the — mixing thus affects the Z-LZP coupling; the Higgs 
mass which is relevant mainly when M^zp ~ ^h/2 and Mkk which sets the mass of the 
new gauge bosons, in particular Z' and Xg. 

We have in addition, r, parametrizing the amount of bulk GUT breaking, and A the 
cut-off scale of the effective 4D theory. Both enter the expression of the X — Xs mixing 
(see Eq. [2T|) thus the LZP coupling to tiXs- The mass of W should be important when it 
can contribute to a coannihilation process near a resonance. Finally, masses of KK fermions 
should be mostly relevant when they are near M^zp and contribute to coannihilation. 

All the formulae for the calculation of masses and couplings needed for the implementa- 
tion of the model into CalcHEP [12\ are listed in the appendix. 

The Z coupling to the LZP is induced via Z — Z' mixing as well as via the mixing with 
9'^, the LH KK neutrino that belongs to the 5 = 1/3 multiplet that contains the SM LH top. 
The Z' coupling to the LZP is of order 1 and proportional to Qiq. The Z — Z' mixing therefore 
increases with qiq and the resulting induced LZP-Z coupling goes as gf^. The — mixing 
goes approximately like Mlzp/M^, thus is large when the the left-handed neutrino is light. 
This corresponds to cc^ = —0.1. Furthermore, the mixing increases significantly when the 
mass of the LZP approaches that of the z/^. Figure [T71 shows how the two components of the 
Z-LZP coupling vary with Mlzp- The component which arises from ul — vr mixing clearly 
dominates, and steeply increases when Mlzp > 1 TeV. 

The Higgs coupling to the LZP is also suppressed by a KK mass. Explicitly, in the limit 



where M^y'^ ^ Mkk and the term under the square root is there only if c,y^ < —1/2. Exact 
expressions used in the code can be found in the Appendix. This coupling increases with 
Mlzp so the largest contribution of the Higgs exchange is expected for a heavy LZP. On the 
other hand, the Higgs contribution is important for Mlzp ~ Mh/2. 

5 Relic density 

In this section, we compute the LZP relic density and explore the parameter space of the 
model. The new features of our calculation compared to Ref. [HI E] are the following: 

• We include all annihilation and coannihilation channels involving level 1 KK fermions 
of the third generation (with the exception of the KK partner of the right-handed top 



M,/ » Mj 



LZP^ 




(12) 
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Figure 17: LZP-Z coupling as a function of the LZP mass. We show separately the contri- 
bution from Z-Z' mixing (dash-dot) and vr-v'i^ mixing (dash) as well as the total coupling 
g'l (full), a) Mkk = 3 TeV, ^lo = 0.3; b) Mkk = 6 TeV, ^lo = 0.3 (blue/dark) and 
^^lo = 1.2 (green/light grey). Lower (upper) curves are for c^/^ = 0.4 (—0.1). 

that is heavy) that belong to the multiplet. We include as well the exchange of level 
1 KK gauge bosons, W, Z' and vector leptoquarks, Xg- 

• We include all possible resonances in either annihilation and coannihilation channels 
and perform precise evaluation of annihilation cross-section near resonances, in partic- 
ular the Higgs resonance which was neglected in [HI [7] . 

• We exactly solve the Boltzmann equation for the LZP number density. Specifically, 
we do not rely on the non-relativistic approximation a = a + hv^ (which fails near the 
resonances) . 

We achieve this using micrOMEGAs 2 . 0. We have rewritten the model in the CalcHEP no- 
tation, specifying the new particles and their interactions. All details can be found in the 
appendix. Since we have a Z3 symmetry rather than a R-parity, special care had to be taken 
to specify the particles that could coannihilate with the LZP. We included all level one KK 
fermions. The diagrams that contribute to each (co)-annihilation process are chosen auto- 
matically according to the interactions specified in the model file. We have also written a 
module for the direct detection cross-section. This only includes the dominant contribution 
that arises from the Z-exchange diagram. The value of gu-, Eq- ( 0^21) . in the LZP model is 
indeed typically too small for the Higgs exchange to contribute, according to the analysis of 
Part I. 

Unless otherwise noted, we consider the following range of values for the free parameters 
of the model: 
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3 TeV< Mkk < 10 TeV, 115 < < 500 GeV, 2Mkk < A < QMkk 

-0.1 < cp^ < 0.9, < r < 0.2, 0.3 < ^lo < 1.2 

The central value Cp^ = 0.4 is fixed by the zero mode of the LH top quark cl and we allow 
deviations from bulk breaking effects (±0.5). The lower value for M^k is constrained by 
EW precision tests and direct detection experiments. The upper value is set arbitrarily to 
concentrate on models that are relevant to low energy phenomenology. 

A light LZP could contribute to the invisible width of the Z. We have taken this con- 
straint into account, V^™ < 1.5 MeV, and found that it plays a role only for the near maximal 
value of (7io ~ 1.2, and Mkk ~ 3 TeV. 

5.1 LZP annihilation 

We first look at self-annihilation before analysing the impact of various fermion coannihila- 
tion channels. We consider the range —0.7 < Cy^^ < which leads to a LZP mass ranging 




MLZp(GeV) MLZp(GeV) MLZp(GeV) 



Figure 18: Relic density of the LZP without coannihilation. The blue full curves correspond 
to (710 = 0.3 and the red dashed curves to Qiq = 1.2. Lower (upper) curves are for c,>^ = 
0.9(— 0.1). The WMAP preferred region corresponds to the horizontal lines. Here, the mass 
of all other KK fermions is 0.9 x Mkk- 

from 1 (2.4) GeV to 2 (4) TeV, for Mkk=^ (6)TeV and A = 2Mkk- Figure M shows 
the behaviour of Qh"^ as a function of the LZP mass for two extreme values of the 5*0(10) 
coupling, gio = 0.3, 1.2 and Mkk= 3,6, 10 TeV. We see respectively the effects of the Z, 
the Higgs and Z' resonances. At large LZP mass, the main annihilation channels are into 
W pairs, heavy fermions (6, t) through Z' exchange and into top quarks through t-channel 
exchange of Xs- 
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Figure 19: Same as Fig. [18] with a) Mh = 115 GeV, b) Mh = 500 GeV. 
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Figure 20: Relic density of the LZP for Mh = 500 GeV, gio = 1.2, c^^ = 0.1, r = 0.1, A = 
2Mkkj corresponding to 10"^ < Qz < 7.10"^ and 6x = 3.2 10""^, with (full) and without 
(dash) the contribution of Xs- 

5.1.1 Annihilation into top 

Models with Pati-Salam SU{4:) x SU{2) x SU{2) gauge group instead of S'O(IO) share most 
of the properties that we have discussed. As far as dark matter annihilation is concerned, 
the main difference comes from the absence of the Xs gauge boson. For this reason, we 
estimate separately the contribution from the t-channel diagram with Xs exchange. It is 
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particularly important near the threshold where the s-channel contribution of neutral gauge 
bosons is not so large. In Fig. [20l we compare Qizph^ with and without the contribution 
from Xs- The difference can reach an order of magnitude. As one moves in the TeV region, 
the shift is more modest as the annihilation into W pairs becomes much more important. 
The net impact of ignoring the contribution of Xs would be to increase the lower bound on 
Mizp from 250 to 600 GeV for this choice of parameters. 

5.2 Coannihilation with KK fermions 




Figure 21: LZP relic density for Mkk = 6 TeV, gio = 0.78, c^'^ = 0.4, r = 0.1, A = 2Mkk, 
corresponding to 5.10"^ < gz < 10"^ and 6x = 3.10"^, including the contribution of coan- 
nihilation channels. The NLZP is a) and b) {tl^vl}- The NLZP-LZP mass difference is 
fixed to A = 0.3 (black/full), A = 0.1 (red/dash) and A = 0.05 (blue/dash-dot). 



As we said, all KK fermions belonging to the multiplet containing t/j, except the KK 
mode of t/j, are expected to be light (compared to the KK gauge boson mass) and close to 
the LZP mass. Among all possible coannihilation channels, we first examine those with KK 
leptons. The couplings of KK leptons r^, r^, z/^ to the Z are much enhanced compared to the 
LZP which is the only one to be suppressed by mixing angles. Therefore, we expect a large 
impact of coannihilation channels for Mlzp ~ 100 GeV where the Z exchange dominates. 
On the other hand, the coupling of KK leptons to the Z' are only slightly suppressed as 
compared to the Z'-LZP coupling and coannihilation effects are more modest when the Z' 
exchange dominates as illustrated in Fig. [211 The coannihilation channels are dominated by 
the NLZP pair annihilation into fermion pairs. Fig. [22] shows the impact of coannihilation for 
two sets of parameters where Qizph"^ is much above WMAP if one considers only annihilation 
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Figure 22: LZP relic density versus the NLZP-LZP mass difference, A, when the NLZP is 
(red), tr (black). Here, gio = 0.78, Mh = 500 GeV c^^ = -0.55, Cj^'^ = 0.4 and 
C/ = 0.1 for all other KK fermions. 

processes. The mass of the LZP is fixed to 200 GeV. We vary the mass of either the left- 
handed KK leptons or the right-handed lepton, tr. For left-handed leptons, the relic density 
drops within the WMAP range for A^, ~ 4% while a much smaller mass splitting is needed 
for TR {Ar = 1%). 

We have also considered coannihilation channels with KK quarks, specifically 6^,^^ since 
6ij has no Z3 charge. Coannihilation channels are dominated by the NLZP pair annihilation 
into fermions and gauge bosons, in particular channels involving quarks and gluons. Because 
of the strong coupling involved, one does not need such a small mass splitting with the LZP 
as in the KK lepton case. For the parameters of Fig. [22l a mass splitting Ag = 12% is 
sufficient. We ignore coannihilation with KK gauge bosons for several reasons. First, the 
motivation for studying the LZP model is that there is a strong rationale for the LZP to be 
significantly lighter than KK gauge bosons. Second, as we have seen above, LZP's in the 
TeV range have already a very low relic density, coannihilation would further reduce this 
relic density. Third, we made the symplifying assumption that all KK gauge bosons are 
degenerate, this approximation is not appropriate in the case of coannihilation where the 
mass difference with the LZP should be known precisely. 

5.3 Direct detection and Summary 

As pointed out in P,[7], the LZP models have large rates in direct detection. In Figure 1231 
we plot cfn-LZP as a function of the LZP mass for the same range of parameters as in Fig[T8| 
compared to the CDMS and XENON limits [131 [H]. For Mkk = 3 TeV, the Z-Z' mixing is 
large and the LZP-nucleon cross-section is above the XENON limit for most of the parameter 
space. From these figures, we conclude that it might be difficult to satisfy both the WMAP 
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Figure 23: LZP scattering on neutron for the same parameter values as Fig. [THJ Solid blue 
and red dashed curves correspond respectively to gio = 0.3 and gio = 1.2. Lower (upper) 
curves are for c,>' =0.9 (-0.1). The two dash-dot lines are the CDMS and Xenon constraints 
[I31[T1]. 

upper bound on the relic density and the limit from direct detection. 

To investigate this more carefully, we scan randomly, generating 10^ models within the 
parameter space 3 TeV < Mkk < 6 TeV, .3 < gio < 1.2, —0.1 < c,>^ < 0.9 and —0.7 < c^j^ < 
0. Here we fix A = 2Mkk and r = 0.1. We keep only models for which flLzph"^ < 0.13. The 
predictions for (Ju-lzp range from ~ 10^'^ to ~ 10^^ pb. Most of these values are already 
excluded by XENON, see Fig. [2H in particular for light LZP's where annihilation proceeds 
through Z or Higgs exchange. Only models for which M^zp > 400 GeV survive. The full 
set of parameters that we have scanned all predict large signals for direct detection and the 
whole parameter space, at least for Mlzp < 1 TeV, will be probed within a few years as this 
requires only a factor 6 improvement in the best limit. If we had scanned up to KK mass of 
10 TeV, the lower limit on cTn-LZP would be roughly the same and we would reach similar 
conclusions once we impose the WMAP constraint. 

The presence of coannihilation channels changes the picture since coannihilation helps 
reducing the relic density but does not affect the elastic scattering cross section. Therefore, 
models with weak couplings to the Z that are still allowed by direct detection experiments 
can also satisfy the upper limit on the relic density. For example, M^zp = 400 GeV, Mz' = 6 
TeV and C/>' = 0.4 (0.9) lead respectively to = 3.7 (2.6)10"'' pb, below the Xenon limit. 
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Figure 24: a) Relic density of LZP for a scan over input parameters, 3 < M^k < 6 TeV, 
0.3 < gio < 1.2, -0.1 < c^^ < 0.9, < r < 0.2, Mkk < A < 2Mi^i^. Here q = 0.1 for all 
other KK fermions and Mh = 500 GeV. Green/light grey dots show all models included in 
the scan, red/grey crosses show only models that also satisfy the XENON bound while black 
crosses show models that satisfy in addition the WMAP upper limit. The WMAP preferred 
region corresponds to the horizontal lines, b) Same as above except that the masses of the 
KK leptons are fixed to 1.05Mlzp- 
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7 Appendix 

Here we give the details of the implementation of the LZP model [6l [7] within CalcHEP and 
micrOMEGAs. 



7.1 Particle Content 

We only consider the first level KK fermions that are in the same multiplet as the SM right- 
handed top quark, which has B = —1/3. We ignore all other KK fermions of the third 
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generation as well as those of the first and second generations. Explicitly we include 



ti, bL, bn, tl, ul, tr, ur = LZP 



(13) 



where we use the shorthand notation = t^P ... The first KK level of Ir is not included. Ir 
has different boundary conditions to provide the zero mode of the RH top quark and its first 
KK mode is in the multi-TeV range. We include of course the zero mode of the right-handed 
top, t^'*. Although we do not include z>^ of the 5 = 1/3 multiplet, this particle plays a role 
because it mixes with the LZP via the top Yukawa coupling. This influences the coupling 
of the LZP to the Z and to the Higgs. We will use the c-parameter associated with it, c^'^. 
The 5D Yukawa coupling of the top is related to the c-parameters of the zero mode of the 
top quark, 

kX, = ^==^==== (14) 



We choose cl = 0.4, cr 



^{1-2cl){1-2cr) 
-0.5 leading to kX^ = 1.58. 



Table 1: KK states included in the model 



name 


Charge 


Z3 


name 


Charge 


^3 


name 


Charge 




tl 


2/3 


-2/3 




-1 


-1/3 


Xs 


2/3 


-1/3 


bL 


-1/3 


-2/3 







-1/3 


Z' 








tR (SM zero mode) 


2/3 





Tr 


-1 


1/3 


w 


1 





bR 


-1/3 










1/3 









The new gauge bosons of SU{2)l x SU{2)r x U{1) are Z' and W (with dominantly 
right-handed couplings to fermions). The two neutral gauge bosons mass eigenstates are a 
mixture of the standard model Z and the first KK mode of the Z' (Z') 



(15) 
(16) 



Z = cos6zz'ZsM + si'^dzz'Z' 
Z' = -sinOzz'ZsM + cosOzz'Z' 

Since the lighter state has to be dominantly the standard model Z, sinOzz' ^ 1, and we will 
ignore the contribution of terms that are suppressed by mixing except when the leading term 
is absent. For charged gauge bosons, the lighter state is denoted by W and is dominantly 
the standard Wl while the heavier state W is dominantly Wr. 

W = cos Oww'Wl + sin 9 ww'Wr 
W' = — sin 9 WW + cos 9ww' Wr 



(17) 
(18) 



Again we ignore the contribution of terms suppressed by mixing, sm.9ww' << 1- The mixing 
angle depends on the overlap between the wave function of the Higgs and that of Z' or W 
denoted ^f, (see Eq. El) 

(19) 
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e 



ZZ' 



910 '^SwCw r>H 

-yz' 



, (20) 
'10 e \Mz'' ^ ' 

In addition, we include tlie leptoquark gauge boson of 5*0(10), X^, which is coloured and 
has electric charge 2/3. The mixing angle Ox between the two 5*0(10) gauge bosons X' and 
Xs appears in the tuiXg vertex and reads (see Ref. [6l [7j for more details) 

9x = kX,{v/Ay (21) 

Here r parametrizes the amount of bulk GUT breaking and will be chosen ~ 0.1. 

7.2 Parameters 

We list the parameters of the model in Table O Unless otherwise noted we take Mxg = 
Mz' = M\Yi = Mkk- The masses M(c, M^k) in TeV of the KK fermions are approximated 



Table 2: Parameters of the model 



Parameter 




Parameter 




Mkk 


Mass of KK gauge bosons 


Cr 


c-parameter of the RH top quark 




\og{Mpi/MKK) ~ log(10i5) 


cl 


c-parameter of the LH top quark 


r 


SO (10) bulk breaking 




c-parameter of the LZP 


A 


Cut-off of RS theory 




c-parameter of r/j 


9io 


4D SO(IO) coupling 




c-parameter of and 








c-parameter of and Bl 


Mh 


Higgs mass 




c-parameter of Br 








c-parameter of 



by (see Fig. 1 of Ref. [7]), 

M(c, Mkk) x 2.405 



M 



KK 



X 2^/a{c) (a(c) + 1)6-'^'"^'^ 



ifc<---. 



Mkk X 2A/c+i if - ^ 



TT 



Mkk X 2(1 + ^) 



if c > 



- e < c < --- 
~ 4 

1 
4 



(22) 



where a(c) = |c+ ||. The numerical value 2.405 arises from the solution to the eigen value 
problem in the Randall- Sundrum geometry where the wave functions are given by Bessel 
functions. 

5D fermions lead to two chiral fermions in 4D. The two chiralities have different boundary 
conditions on the TeV and Planck branes. The LH helicity of the LZP turns out to be 
localized near the Planck brane while the RH helicity is near the TeV brane (see Ref. [7J for 
more details). Therefore, only the RH helicity couples significantly to the KK gauge bosons. 
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7.3 Wave functions 

7.3.1 Fermions: SM tR and KK ur 



rc\2-c 



e^=(l - 2c) 



1/2 



fLZp{z,c,m) = {ze'^Y 



5 e 



erc(i-2c) _ I 
Ja{c){me~'''') 



where 6(c, m) 



N'^{c,m) cK-i/2 

c>-l/2 



7.3.2 KK gauge bosons 



2r. 



{Ja(c){m) + 6(c,m)F„(c)(m))^ 
(^a(c)±i("^e~''=) + 6(c, m)y„(c)±i(me~'''=))^] 



hga{x) = 



='-c/2. 



Ji(xe-''=) 
yi(xe-^=) 



7.3.3 Higgs boson 



2rc 



7.3.4 Wave function overlaps 



^{(C1,C2) 




/ {ze^-^y^ fga{z, 2A05)ft^{z, C2)fLZp{z, ci, 2.405) 


52(Cl,C2) 




/ (^e'-^)-' fga{z, 2A05)fLzp{z, c,, 2.405)/izp(z, C2, 2.405) 






[ dz {ze^'^y^ fgaiz, 2A05)flzp{z, c, 2.405) 






f dz{ze'-^y'Ua{z.2AQh)fl{z,c) 


yz' 


= a/ rcC- 


[ dz{ze^^y'f,,{z,2A05)f^{z) 
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Table 3: Couplings of KK fermions to neutral SM gauge bosons 



Vertex 

hbLZu, 
bRbnG^ 
bRbnAn 
bRbnZ^ 
bnbZf, 



Value 

9s7^^ 

e 



(2 l.^w) 



CwSw 



( 2 + l^w) 7m 



3cw 



2swCW •' '^^ 



/c«M^7.(l-75) 



Vertex 

TrTrZ^, 

tlTlA^ 
tlTlZ^ 
t^t^lZ^, 

^RVrZ^i 



Value 

cw 

e 

i^z"7M(l + 75) 



Table 4: Couplings of KK fermions to W 



± 



Vertex 


Value 




b~RtW^ 

htiW- 


\{gR,eww'i,{i + 75) + ^^^7m(i - 75)) 

\9RiOww'lii{'^ + 75) 


9r^ = 9w\m9if{cbR-, cr) 

9Ri = 9lQ^/^92fic^.R, Cth) 



7.4 Neutrino mixing and coupling to the Higgs 

The large Yukawa coupling between the 16 with S = 1/3 containing t° and the 16 with 
B — —1/3 containing t*^ generates a mixing mass term between i/^ and the LZP. This induces 
LZP couplings to the Higgs and the Z after diagonalization of the mass matrix 



LZP 'n^URu'j^ 







(35) 



where Mlzp = M{cur, Mkk), m^'^ = M{cy>^, Mkk), m^Ru'^ = mtfcLfcnivn- 
Cf^l if c/>-l/2 and Cf^^/fcf if c/ < -1/2 and /c = 



1 - 2c 



(36) 
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Table 5: Couplings of SM particles and KK fermions to new gauge bosons 



Vertex 
IvrXs 

VrUrZ'^ 

ttzi 

bhZ'^ 
HZ.Zl 

bRbuZ'^ 
W^W^Z'^ 

btW'~ 



w^m^z 



Value 

1^i7m(1 + 75) 
f^x7M(l-75) 
^7.(1 + 75) 
^7^(1 + 75) 

^7^(1+75) + 1^7m(1-75) 
1^7^(1-75) 

Sw^w 

^7/.(l - 75) 



^7/.(l - 75) 
%^7/.(l - 75) 
%^7/.(l - 75) 

#7.(1+75) 



V2i 



Vtb{—Oww')lii 



^7m(1 + 75 



sw 



9ri 



-Oww'hfii^ -75) 



f^7/.(l+75) 



e 



V2! 



Sw 



75 j 



91 = 9io\/^9if{cuR,CR) 
9i = 9iG\/^9if{Cu'^,CvR) 
9z' = 9w\f^9z'icun) 

9t = 9W\/^9z'i'^rR) 
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9tR = 9w^9z'icR) 
9tL = -9w\/^9z'{cl) 
9hz = 



9i 



9^1 — -9io^/^9z'{ctL] 



9z' — ~'^9w\/^9z'i'^b 



from Wl — Wr mixing 
from Wl — Wr mixing 



The LH and RH helicities of the LZP (respectively the mass eigenstates of MlMy and M^Ml) 
are 

(z/i)l = cos 6'lz>/j + sin (37) 
{vi)r = cos 9 rVr + sin Orv'^ (38) 

where vr and z/^ are the 5D KK partners of vr and i/^, which have respectively LH and RH 
chiralities (see Ref. [7j for more details). Therefore, the LZP coupling to the Higgs is given 
by: 

9h = fcL fcR sin Ol cos Or^^^^ (39) 
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where we used 2kXr, = fc^fcR- All fermions whose coupling to the Higgs are also induced 
via mixing with a heavier KK state (tr, tl, i'l, ti, bi) have a coupling gujf given by Eq. [39l 
However, in the mixing angle one should substitute the following masses for the fermion / 

mi = MxKicf, Mxs) rn2 = Mkk{cu>^, Mxs) = ^tfcLfcR^f (40) 

Here we assume that all the heavy states which mix with level one KK fermion have a 
common mass equal to Mxi^(c,>^, Mxg)- The mixing angles are 

/ (-m? + - + yA)A ^ ^ 

{-ml + mj - mj + y^) / (-mf + - mj + v^)^ ^ ^'"^ 
cosUr = i H — [4:1) 

2m2m3 \ 47712771^ I 

where mi = Mlzp, ^2 = m^' , rris = m^^^' and A = —Amlml + (mf + m^ + m|)^. 



7.5 The Z-LZP coupling 

The Z-LZP coupling ((7^) has one component due to Z — Z' mixing, g^^i, as well as another 
component from i^'iiyR mixing, g^^^^- Both components are comparable although small. 
Explicitly, 

yz — Hzz' ^ yuLUR 

9Tz' = 9z' dzz' 

gll,^ = g^z sm' 9r (43) 

where g"^ = e/2swCw is the Z coupling to z/£ and sinOiiR is defined in EqJHl For rriy^ <^ 
m,^^, sin "C svhOr and we can neglect the coupling of the LH helicity of the LZP. 9zz' is 
given by Eq. [20] and gz' is the Z' coupling to the LZP, (see Table E]), 



St* 

9z' = 9w\i ^9z'{cur) (44) 
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